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Evolution from a Bose-Einstein condensate to a Tonks-Girardeau gas: 

An exact diagonalization study 
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We study ground state properties of spinless, quasi one-dimensional bosons which are confined 
in a harmonic trap and interact via repulsive delta-potentials. We use the exact diagonalization 
method to analyze the pair correlation function, as well as the density, the momentum distribution, 
different contributions to the energy and the population of single-particle orbitals in the whole 
interaction regime. In particular, we are able to trace the fascinating transition from bosonic to 
fermi-like behavior in characteristic features of the momentum distribution which is accessible to 
experiments. Our calculations yield quantitative measures for the interaction strength limiting the 
mean-field regime on one side and the Tonks-Girardeau regime on the other side of an intermediate 
regime. 
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One-dimensional (ID) delta interacting bosons reveal 
remarkable similarities with non-interacting fermions 
when the interaction between the particles is strong [l|. 
The tremendous experimental progress in the field of cold 
atoms has recently allowed for the realization of such 
strongly interacting ID bosonic systems @, S fl Q ■ The 
opposite regime of weak interactions is well described by 
the Gross-Pitaevskii theory [Q, |7|. Besides the ground 
state properties, the dynamical behavior |8| in these two 
limiting regimes is very different and the excitations fol- 
low the Luttinger liquid theory [3, 1 1011. Moreover, an 
intermediate regime is distinguished [ll| by its charac- 
teristic phase fluctuations indicating the onset of corre- 
lations between the bosons. In the homogeneous ther- 
modynamic limit an exact solution covering all regimes 
has been introduced by E. H. Lieb and W. Liniger [12j |. 
which is the basis of many contemporary approaches 
@) LUa, LLJ, [Da, UM LLZl ■ Correspondingly, most theoretical 
studies assume large particle numbers. However, gases 
with strong correlations could only be realized experi- 
mentally with a small number of particles so far. 

In this article we concentrate on small systems with few 
particles where the interaction strength between the par- 
ticles can be tuned by the transverse confinement [5J. In 
particular, we study the influence of interactions via the 
pair correlation function which clearly indicates the lim- 
its of the mean-field (MF) or Bose-Einstein condensate 
(BEC) regime, an intermediate regime and the Tonks- 
Girardeau (TG) regime. The discrimination of these in- 
teraction regimes is an important question with relevance 
for current experiments. We show that the transitions 
between the BEC, the intermediate and the TG regime 
can also clearly be distinguished in the evolution of the 
momentum distribution. The detailed analysis of the cor- 
relation function and of the momentum distribution is a 
central method in the whole field of quantum gases. 

We consider spinless bosons (e.g. 87 Rb atoms 



with frozen spin degrees of freedom) confined in 
a three-dimensional cigar shaped harmonic trap 
V ex t(r) = \mJ 1 _ (x 2 + y 2 ) + \mw 2 z 2 \ m is the mass of 
the bosons, and lo z and u± are the axial and transverse 
angular frequencies. The transverse confinement is much 
stronger than the axial confinement lu z <C oj±. The 
bosons are assumed to interact via a delta potential 
MntG^ - f'l) = a? S(r — f'), where a s is the s-wave 



scattering length. 

The system becomes quasi one-dimensional when the 
transverse level spacing Hu>± is much larger than the axial 
level spacing ftw z and the three-dimensional interaction 
strength Usd — 47r f ,° 3 (l z , I A.'- oscillator lengths of the 



axial and transverse direction, L = «/—*=-). Under these 
conditions the transverse motion in the ground state is 
restricted to zero-point oscillations. Therefore, the many- 
particle Hamiltonian reads 



H = hw z ^2 



2 ' " ; ' / ' 



\u m 



22^jkia\a]aia k , 



ijkl 



(i) 

where a\ (ctj) is the bosonic creation (annihilation) 
operator for a particle in an energy eigenstate (fn of 
the axial harmonic oscillator. We have neglected the 
constant zero mode energy of the transverse oscillator 
potential, hjkl are the dimensionless interaction in- 
tegrals I i jki=l z J™ 00 dz(t>i(z)(t)j(z)<}>k(z)<fii(z). The es- 
sential parameter to characterize the system is the 
one-dimensional interaction strength [/id- It results 
mainly from an integration over the transverse direc- 
tions. In addition, the strong vertical confinement leads 
to a modification of the s-wave scattering length [18j |: 
a c fr = a s /(l — lA6a s /v2lx). In this paper we restrict 
ourselves to confining frequencies relevant to the exper- 
iments by Kinoshita et al. [5(, co± — . . . 2tt x 70.7kHz, 
resulting in corrections to a s oia s < a e e < 1.16a s . These 
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FIG. 1: Particle density of five bosons for different interac- 
tion strengths [/id along the axial direction (x = y — 0). The 
density becomes flatter and broader with Uiu- In the strong 
interaction regime it developes oscillations. At Uio — 20hu) z 
the density of the bosons perfectly agrees with the density of 
non-interacting fermions. 



FIG. 2: Pair correlation function of five bosons for different in- 
teraction strengths [/id along the axial direction (x — y — 0). 
One particle is fixed at position z — 0. The distribution flat- 
tens and forms a hole at coinciding particle positions, z — z'. 
The correlation function clearly indicates the transition be- 
tween the three regimes (see text). 



values of the transverse frequency are far from the con- 
finement induced resonance discussed in Ref . [181 . Then, 



U 



ID 



U3D g cff 

2-k a s 



2y/mhuj z a c f[UJ± , indicating that the ef- 



fective interparticle interaction can be tuned by the 
transverse confinement. The Hamiltonian {1} is diago- 
nalized in the subspace of the energetically lowest eigen- 
states of the non-interacting many-particle system, con- 
sisting typically of up to 150000 basis states. In the fol- 
lowing we will concentrate on results achieved for five 
bosons. 

We start our discussion with the particle den- 
sity p(z) = (# t (z)*(z)) (#(z): field operator) which 
is shown in Fig. [T] At small interaction strength 
the density reflects the conventional mean-field be- 
havior and p(z) m N4>mf(z) 2 . In this regime all the 
bosons condense into the same single-particle wavefunc- 
tion, ^>b{z\, ■■■,Zn) ~ ris=i ^MF^iji and thus the many- 
boson system is well described by <^mf(z), which solves 
the Gross-Pitaevskii equation |6|,|7J]. The system reacts to 
an increasing repulsive interaction with a density which 
becomes broader and flatter 



ia 17, ll 



{j. In the 

strong interaction regime density oscillations appear (see 
e.g. the curve at Um = 8Huj z in Fig. Q]) and with fur- 
ther increasing Um the density of the bosons converges 
towards the density of five non-interacting fermions, 
Pf(z) — Y^i=o ^?( z ); as predicted by Girardeau [l(. Both 
densities are in perfect agreement at Um — 20hcu z indi- 
cating that the limit of infinite interaction is practically 
reached. Thus, the density oscillations reflect the struc- 
ture of the occupied orbitals in the harmonic trap. In 
contrast to Ref. 2l| which predicts the oscillations to 
appear one after the other, when the repulsion between 
the bosons becomes stronger, we observe a simultaneous 
formation of five density maxima. These density oscilla- 



tions are absent in mean-field calculations 13|, [l9| . 

Additional insight into the evolution of the 
system with increasing interaction strength can 
be obtained from the pair correlation function 
g i - 2 \z,z') = ffi(z)&(z')$l(z')$t(z)) which is depicted 
in Fig. [3 for different Um- In the weak interaction 
regime where the mean-field approximation is valid 
the correlation function resembles the particle density 
and 5 ( 2 >(z,z')« N{N - 1)4>mf(z) 2 4>mf{z') 2 - The ap- 
pearance of a minimum at [/id = 0.5fiw z marks first 
deviations from this mean-field behavior. The interpar- 
ticle interaction leads to a reduced probability of finding 
two bosons at the same position. These correlations are 
characteristic for the intermediate regime. With increas- 
ing interaction strength the correlation hole increases 
and already at Um = 3huj z the correlation function at- 
tains a form which is typical for a Tonks-Girardeau gas: 
Flat long-range shoulders indicate the incompressibility 
of a Fermi gas. This interaction strength thus marks 
the transition from the intermediate to the TG regime. 
By contrast, the density still exhibits a mean- field 
shape. The correlation function reaches its limiting 
form corresponding to the one of five non-interacting 
fermions at Um — 20hut z . Besides the central correlation 
hole maxima appear which indicate the position of the 
other particles. In our finite size system this limiting 
behavior is reached at smaller interaction strength than 
in homogeneous systems [22| |. Thus, the pair correlation 
function clearly indicates the limit of the mean-field 
regime at small interaction strength (Um = 0.5hoj z ) 
and the transition from the intermediate towards the 
Tonks-Girardeau regime at larger interaction strength 
([/id ^ 3hui z ). Within our calculations these limits are 
not sensitive to the number of particles, which we have 
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FIG. 3: Evolution of various contributions to the total energy, 
Etot.z, of five bosons with increasing interaction strength Uid- 
The energies evolve towards the accordant energies of non- 
interacting fermions [231 ] . The minimum of the kinetic energy 
coincides with the onset of significant correlations and there- 
fore marks the limit of the MF regime, 
decreases in the whole interaction regime, see inset. 



By contrast E^^ z 



FIG. 4: Momentum distribution of five bosons for different 
interaction strengths Uid (j?x = Py — 0). The central peak 
reaches its maximum height at Uid — 3hui z , see inset. At 
[/id = 20^2 the momentum distribution has developed the 
typical form of a harmonically trapped TG gas [2_| which is 
completely different from the distribution of non-interacting 
fermions. 



checked for up to N = 7. 

We note, that due to the singular shape of the inter- 
action potential the correlation function exhibits kinks 
at coinciding particle positions, z — z', 22|, |24| which 
are not resolved in Fig. [__ |25fl. In recent experiments 
the correlation function of three-dimensional ultracold 
atomic systems has been obtained by analyzing the shot 
noise in absorption images [2__, [27J. Its value at z = z' 
1J, |l__, LL7J determines, e.g., photoassociation rates [28[ 
and the interaction energy. The latter is given by 
E int = ^p / dzg( 2 \z, z) and is depicted in Fig. ___ Its 
behavior is similar to the homogeneous system due to 
the short range of the interaction potential [1__] . Never- 
theless, measurements of the whole pair correlation func- 
tion are tedious. A quantity which is easier accessible to 
experiments is the momentum distribution. In the fol- 
lowing we demonstrate that the transition between the 
different regimes discussed above can also be obtained 
from this quantity. 



We first discuss the kinetic 



TP - Ml 

-C'kin.z — o™ 



energy, _*,___._. - 2m 
which is proportional to the width of the momentum 
distribution. The limit of the mean-field regime at 
[/"id ~ 0.5/kj z is clearly visible in the onset of a min- 
imum of the correlation function (Fig. [2]). At the 
same time the kinetic energy changes dramatically, see 
Fig. [3J In the weak interaction regime it can be approx- 

frM F(z) The 



particle distances lead to an increase of the kinetic 
energy which can be read from the exact expres- 



<)zl 



&b(zi ■■ -zn) 



By 



imated by £ kin ,z « E^ z = Nfa f dz 
flattening of the density therefore results in the ini- 
tial decrease of the kinetic energy. However, this ef- 
fect is in competition with the development of short 
range correlations in the intermediate interaction regime. 
Strong variations of the wavefunction at small inter- are responsible for the increase of the kinetic energy, i.e. 



sion .Ekin.z = N^J dz x . .. dz N 
contrast, the mean-field kinetic energy, which is only sen- 
sitive to variations of the density, decreases in the whole 
interaction regime, see inset of Fig. [__ Therefore, the 
minimum of the exact kinetic energy clearly marks the 
limit of the mean-field regime and the dominance of in- 
terparticle correlations. With increasing particle number 
N the minimum of the kinetic energy slightly shifts to- 
ward smaller values of [/id with a limit at Uyd — 0-5Huj z 
at large N. 

We mention that the potential energy (Fig. [3]) is 
proportional to the width of the particle density, 
■Epot.z = ^mojl(z 2 ). In the strong interaction regime 
both quantities reach its limiting value of the non- 
interacting fermionic system. In experiments [5( this has 
been used as an indication for the Tonks-Girardeau limit. 

While the width of the momentum distribution in- 
dicates the limit of the mean-field regime the evolu- 
tion of its central peak marks the transition towards 
the Tonks-Girardeau regime (Fig. [4|. This central re- 
sult of our studies does not depend on the number of 
particles. The momentum distribution is defined as 
p(pz) = (IF (p z )Il(p z )) where II (p z ) annihilates a boson 
with momentum p z . Corresponding to the flattening 
of the density, the height of the central maximum of 
the momentum distribution increases at small interac- 
tion strength. Already in this regime, high momentum 
tails develop due to the kinks in the wavefunction at co- 
inciding particle positions, Zi = Zj 3Qj, __1|- These tails 
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FIG. 5: Occupation number distribution of five bosons for 
different interaction strengths Uyd- With increasing [/id the 
bosons leave the ground state and occupy excited states. 
Single-particle states with even parity are comparatively 
stronger populated than those with odd parity. 



{pi), at small interaction strength above U\t> = 0.5huj z . 
However, the width of the central peak is still shrink- 
ing in this regime, corresponding to a further broaden- 
ing of the density. By contrast, the formation of the 
correlation hole at z = z' leads to a redistribution from 
low towards high momentum states. This effect domi- 
nates above Uuy = Shuj z , when the growth of the density 
width slows down. At this point the height of the cen- 
tral peak has reached its maximum. This coincides with 
the transition behavior visible in the correlation function. 
Therefore, the interaction strength at which the central 
peak of the momentum distribution reaches its maximum 
height marks the transition towards the Tonks-Girardeau 
regime. Within our calculations the value of U\y> — 3fiw z 
marking this transition point is independent of the par- 
ticle number N . The experiments of Tolra et al. [3j 
(Uir, ~ 4.91 hu z ) and Kinoshita et al. [5j (U\t> up to 
« 15.4 hu> z ), therefore, both have been performed in the 
Tonks-Girardeau regime. The height of the central peak 
at its maximum is approx. 30% larger than at small inter- 
action strength ([/id ~ 0) and about 20% larger than at 
large interactions ([/id = 20fiw 2 ). This contrast increases 
with increasing particle number. Moreover we mention 
that the momentum distribution at U\u = 20huj z per- 
fectly agrees with exact results obtained in the limit of 
infinite interaction strength [29j |. 

Finally we discuss the occupation number distribution 
{rii) = (a]ai) of the ID harmonic oscillator states which 
is shown in Fig. \5\ With increasing interaction strength 
[/id the bosons leave the ground state and occupy ex- 
cited states. At [/id = 20huj z the distribution is similar 
to the distribution shown in [32j for [/id = oo. How- 
ever, we observe a stronger population of single-particle 



states with even parity compared to those with odd par- 
ity. This effect is most pronounced in mean-field calcula- 
tions where occupations of odd parity orbitals are absent. 
The comparatively stronger occupation of single-particle 
states with even parity can therefore be interpreted as a 
remnant of the mean-field regime. 

In summary, using the exact diagonalization method 
we studied the interaction-driven evolution of a quasi 
one-dimensional few boson system. Besides the pair cor- 
relation function we identified the momentum distribu- 
tion as a reliable indicator for transitions of the sys- 
tem between three characteristic regimes, the BEC or 
mean-field regime, an intermediate regime with strong 
short range correlations and the Tonks-Girardeau regime. 
From this we were able to quantify the interaction 
strength for the transitions. The width of the momen- 
tum distribution has a minimum when the interaction 
strength is approximately half as large than the axial 
level spacing of the trap ([/id = 0.5huj z ). This behavior 
coincides with the onset of significant correlations and 
therefore marks the transition between the mean-field 
and an intermediate regime. The central peak of the mo- 
mentum distribution reaches its maximum height when 
the interaction strength is approximately three times 
larger than the axial level spacing ([/id = 3fruj z ) and al- 
ready at this point the pair correlation function attains a 
form which is typical for a Tonks-Girardeau gas. The evo- 
lution of the central peak of the momentum distribution 
therefore marks the transition between the intermediate 
and the Tonks-Girardeau regime. These features of the 
momentum distribution allow a reliable discrimination 
between the three regimes. We are aware of the limita- 
tions of our method to small particle numbers, however, 
we want to point out that the method of exact diagonal- 
ization is capable to reveal the basic microscopic mecha- 
nisms of quantum gas systems which often determine the 
physics of larger systems. 

We thank S. Reimann, M. Ogren and H. Monien for 
fruitful discussions. 

Note added in proof. Recentl y, w e became aware of 
related work by S. Zollner et al. 331.134 1. 
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